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2i ■ Abstract 

["t I ■ It is sliown that there exists a range of parameters in which gravitational 

h»^ , collapse with a spherically symmetric massive scalar field can be treated as 

CN I if it were collapsing dust. This implies a criterion for the formation of black 

holes depending on the size and mass of the initial field configuration and the 

^ ■ mass of the scalar field. 
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I. INTRODUCTION 

An important issue in General Relativity is to ascertain whether a given matter distri- 
bution will collapse to form a black hole. To address one aspect of this question we shall 
consider a minimally coupled massive scalar field with spherical symmetry which starts from 
rest. Massive scalar fields are important in many models of the early universe. If any of 
these fields underwent gravitational collapse then the result would be a primordial black 
hole. 

In recent years, gravitational collapse with massless scalar fields and spherical symmetry 
has been analysed in some detail p|-§[. One result to come out of this is that, when the 
initial profile of the scalar field is fixed in terms of a single parameter, there is a critical 
parameter value. Beyond the critical value a black hole will form 0] . 

Much less is known about the spherically symmetric collapse of massive scalar fields. 
In order to obtain an approximate solution of the Einstein equations we introduce a WKB 
approximation for the field in the limit where the mass of the scalar field is large. To leading 
order, the solution behaves in the same way as inhomogeneous dust. This is known already 
for the linearised Einstein equations, but we find that the WKB approximation is equally 
good in the non-linear regime. 

The dust solution, like the Oppenheimer-Snyder solution, always collapses to a black 
hole. By keeping track of the size of the terms that are higher order in the approximation 
we obtain a sufficient condition for the validity of the approximation and therefore for the 
gravitational collapse of the massive scalar field. This condition can be obtained analytically 
when the initial density has a top-hat form, and numerically otherwise. 

Planck units (in which G = c = h = 1) are used throughout. 

II. THE MODEL 

The most general spherically symmetric spacetime metric can be written in the form 

ds^ = -e^^dr^ + e-2^dr2 + i^M(]^ (1) 

where \E' = \E'(r, r), A = A(r, r), R = R{T,r) and dfi^ = d6^ + sinOdcj)^ is the metric of a 
unit 2-sphere. It is convenient, for our purposes, to rescale the coordinate r to the proper 
time t of a co moving observer at radial coordinate r, 

t= re*(^''^)dr', (2) 

Jo 

hence one gets the metric in Gaussian coordinates, 

ds^ = -dt^ + e-''^dr^ + R''dn^. (3) 
The Einstein tensor components are 

Gu = R'^ [-Re^^{2R!h' + 2R" + R-^R!^) - 2RAR + 1 + R^] (4) 

Grt = -2R-\R! + R!K) (5) 

'e~^^{2RR + R^ + l)-R!^] (6) 



yj 7->7^ J. L 



-2 



' rr 



Gee = sin-2 OG^^ = R[RA + R'A'e^^ + R"e^^ -R + AR- A^R] (7) 



where dt = ■, dr = ' . 

For the matter content, we will take a real minimally coupled scalar field of mass /i, 
whose equation of motion is 

(-V2 + ^2)0 = 0. (8) 

In the spherically symmetric metric, 

- e2^0" + li^ + (A + 2R-^r) - e^^ (A' - 2R-^B!) 0' = (9) 

The stress-energy tensor of the scalar field has components 

Tu = \^' + ^e^V'^ + ^/i^0^ (10) 

Trt = 00' (11) 

T; = \^' + \e'^r - ^/iV^ (12) 

Tee = sin-2 OT^^ = \R\^^ - e^'^^'' - i?<i?) (13) 

It proves convenient to introduce the functions 

k{t,r) = l~e^^R!^ (14) 

m{t,r) = \R{R^ + k), (15) 

so that two Einstein equations can be recast in terms of the first derivatives of these two 
functions, 

k = 87rRR'T[ (16) 

m = A-kR^R'TI - AttR'^RT^ (17) 

These two equations, together with the scalar equation @, form a complete set. A third 
Einstein equation provides a constraint 

m' = inR'^R'Tu - Ar^R^RTrt (18) 

which is important for relating m{r, 0) to the initial data. 

The equations have an approximate solution when the compton wavelength of the field 
/x~^ is much smaller than the radius A of the spherical region where the field is non- vanishing, 

A/i > 1 (19) 

The scalar field has a wave-like solution with nearly constant amplitude, 

0(t, r) = /i"^<l>(t, r) cos(/it). (20) 

The stress-energy tensor becomes 



Ttt = 1$^ - l/i"^** sin 2/it + |/i"^(<i'^ + e^H'^){l + cos 2/it) (21) 

Trt = -l/i"^**' sin 2/it + ^fi-^<l><t>'{l + cos 2/it) (22) 

T/ = -i$2cos 2/it-i/i-^$<i>sin 2/it + i/i-2(<j)2 + e^^$'2)(l + cos 2/it) (23) 

Tgg = \R^ [-$2 cos 2/it - ^-^•Ixi' sin 2/it + i;U-2(<i.2 _ e4A$'2^)(^^ ^ ^Qg 2/it)] (24) 

The form of the stress-energy tensor suggests a trigonometric expansion for the function 
$ of the form 

oo 2m 

$ = $0 + ^ 5] /i"'" ($^„ COS n^it + <l>^„ sin n/xt) (25) 

m=l n=2 

with similar expansions for m and k. More precisely, we take 

$(t,r) = $oi(t, r) cosyUt + /i"2<l>22 cos 2/it + 0(/i"^) (26) 

m{t^r) = mo(t,r) + fi'^mi2{t,r) sm2nt + jj,~'^m22{t,r) cos2ixt + 0{jj^^) (27) 

A;(t, r) = ko{t, r) + fi~'^k22it, r) cos 2/it + 0(/i"^) (28) 

R{t,r) = i?o(t,r) + /i-2i?22(t,r) cos 2/it + 0(/i-^) (29) 

Terms which play no role have been dropped. 



When the expansions are substituted into eqs. ([TBHTTD the leading order terms give 

mo = (30) 

k = (31) 

Rl = 2moRo^ - ko (32) 

m[, = 2txRIr!^^1 (33) 

The first three equations also arise in the Tolman-Bondi solution representing the collapse of 
spherically symmetric dust. We deduce that to leading order the metric will have Tolman- 
Bondi form. 

The next order gives 

A:22(t,r) = 27ri?o(i?o)-'(l - A:o)$o$o (34) 

mi2{t,r) = 'nRlRo^l (35) 

i?22(t,r) = -|i?o$o> (36) 

We will assume that the WKB approximation holds as long as these correction terms are 
less than 1/2 the leading terms. (The effect of changing this number will be equivalent to 
changing the value of /i). The largest correction comes from the mi2 term, therefore the 
WKB approximation should be valid whilst 

mo 2 

This inequality defines a region of the (t, r) coordinate plane. The region is bounded by a 
curve where the inequality is saturated, 

Rof^'o = firriQRQ. (38) 

There is no reason to extend the leading order approximation beyond this curve. Therefore 
the region of validity of the approximation will only consist of points in whose causal past 
the inequality is always satisfied. 



III. TOLMAN BONDI METRICS 

We have seen that a sufficiently massive scalar field behaves like inhomogeneous dust. 
The corresponding metric can be recovered from the definitions (|1^) and (|15]), with m{t, r) = 
m{r) and k{t,r) = k{r) (dropping the 'O'-indices for simplicity), 

ds^ = -dt^ + ^5^%4dr2 + R\t, r)dn^ (39) 
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The metric has Tolman-Bondi form 03. This is a well known class of solutions that can 
be thought of as a collection of free-falling spherical shells. Assuming R' y^ (then none of 
the shells cross), one can analytically integrate equation (^) parametrically 

Tfl 

t{r],r) = to{r) + --j{T] + smT]) (42) 

k2 

Riv,r) = —cos^ (43) 

where to{r) is an arbitrary function. 

We take time-symmetric initial data with 

R{0,r)=r (44) 

R{0,r)=0 (45) 

These initial conditions fix the radial coordinate r and the functions tQ{r) = 0, k{r) = 
2m{r)/r, 

R(r], r) = r cos^ - (46) 






tiv,r) = d^{v + sin r]) (47) 
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where rj runs from to tt. A particular solution is uniquely specified by m{r) and related 
to the initial density profile Ttt{0,r), 

m{r) = 471 r Tu{0,r')r'^dr' (48) 

The function m{r) is the total mass within the sphere of radius r. If it approaches a constant 
value M at large radius then M is the ADM mass. 

A model for the initial density profile which we make use of later is 

p(r) = Tii(0,r)=poe^(^'/")', (49) 

so that the mass function is simply. 



m(r) = M(l-e-("/^)') (50) 

This mass distribution collapses to form a black hole of radius 2M. The centre collapses 
first, and forms a singularity after a proper-time tt/^/SM. This is shown in figure [^ 

Also shown in figure |I] is the event horizon, obtained by solving the equation for radial 
null geodesies 
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±(—-1] R' ---f(r] + sin r]) R'^ (51) 
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dr 
where 

7(r) = - - 1:11 (52) 

'^ ^ 2 2rr ^ ^ 

with the boundary condition R —>■ 2M as r -^ 00. 

Various derivatives of the radial function that are of use are 

R = -[—) tan^ (53) 

77 X 77 

R' = cos^ 7: + 7:1 iv + sin 77) tan - (54) 



The value of $ given by eq. (^TD will diverge if i?' = 0, from which we obtain a restriction 
on the initial mass distribution: 7 > 0. 

IV. THE OPPENHEIMER SNYDER LIMIT 

Oppenheimer and Snyder considered the collapse of a star with constant density, for 
which 

mir) = l^{j) < ^ < ^ (55) 

[m r>X 

In this case the model can be solved explicitly and we can obtain analytic expressions for 
the domain of validity of the WKB approximation. We are especially interested in those 
cases where the domain of validity covers the whole region outside the event horizon. 

A constant scalar field gives a constant density but the discontinuity at r = A causes the 
WKB approximation to break down. It is possible to change the form of the field close to 
the edge to make the field continuous and still have constant density, 

{$0 r < A - /i"^ 

<l>i(r) \-fi-^<r<\ (56) 

r > A 

where $0 = 0(1) and $; = 0(1). 

The WKB approximation breaks down along the curve given by equation (|38D. Let ri^{r) 
parameterise this curve, then from eqs. (^) and (|54D , 



Therefore ?7*(r) is constant for r < X. In this case t{ri^,r) is also constant. The Tolman- 
Bondi solution is exact for r > A, therefore the WKB approximation breaks down along a 
line segment from r = to r = A with fixed time coordinate t*. If this line segment lies 
inside the event horizon of the Tolman-Bondi metric then the WKB approximation must 
be valid outside of the event horizon. This is shown in figure |^. 

The event horizon is given by a curve //^^(r). Outside of the matter distribution the 
event horizon lies along R = 2M, therefore 

2 riEHJr) 2M 

cos = r > X (58 

2 r - ^ ^ 

Now, the end of the line where the WKB approximation breaks down will lie inside the event 
horizon when 

V* > VehW (59) 

which is when 



tan^>i/— -1 (60) 

2 V 2M ^ ^ 

Now, from equation (|57D we have 

tan^fl+tan^^U^./^- (61) 

2 V 2 / 3 V 2M ^ ^ 

The condition that the WKB approximation breaks down inside the event horizon is therefore 

(62) 

If this inequality is satisfied, then the scalar field configuration with ADM mass M which 
is initially constant for r < X will collapse to form a black hole. The arguments presented 
here say nothing about what happens when the inequality is not satisfied. The inequality 
approaches /iM > 3/2 when A is much larger than the initial Schwarzschild radius. 

V. INHOMOGENEOUS MODELS 

For more complicated density distributions it becomes necessary to test the validity of 
the WKB approximation numerically. We will concentrate here on the exponential density 




distribution (^91) given earlier and compare this to the results for the top-hat distribution of 
the Oppenheimer-Snyder limit. 

The aim, as before, is to locate the parameter ranges for which a black hole forms 
because the WKB approximation holds good everywhere outside of the event horizon. The 



method we use is to find the curve on which the WKB approximation breaks down and find 
conditions which place this curve wholly inside the event horizon. 
Differentiating eq. ( pSD gives an equation for r]^{r), 



7-i(3 + cot^|) 
with 77^,(0) = ?7*, given by eq. (^) and 



2 V* 
cos - 



^=/5'-^7'(^* + sinr^.) (63) 



, , 3 rm' , . 

^'<'-' = 2-2^ («^' 

P(r)=(?^y^ (65) 

This equation can be solved numerically using a simple Runge-Kutta method. An example 
is shown in figure |^. 

When the shape of initial field configuration is fixed, the model is parameterised by 
the scalar field mass /i, the length scale A and the total mass M. The method used here 
only depends on the combinations A/x and M/X. The numerical results for two density 
distributions are shown in figure ^. Above the lines, the WKB approximation is valid 
outside of the event horizon and a black hole will form. This condition is sufficient but not 
necessary. The numerical results for the Openheimer-Snyder limit agree with the analytical 
result to six significant figures on the majority of data points. 



VI. CONCLUSIONS 

We have analysed in some detail the regime in which a spherically symmetric massive 
scalar field can be treated as if it where dust. The result is a criterion for the formation of 
black holes depending on size and mass scales. The criterion depends to some extent on the 
form of the density distribution. 

Outside of the WKB regime the solution of the Einstein equations is a more difficult 
problem. The most studied case is the massless limit, where echoing behaviour was discov- 
ered at the point which separates the cases which form black holes from those which do not. 
This point lies on the horizontal axis in figure |^. It is quite possible that the line which 
bounds black hole formation continues up from this point. It certainly must lie under the 
lines shown in the figure. 

In some infiationary models of the early universe, the end of the infiationary phase is 
followed by a period dominated by the energy of a massive scalar field P-ITll. It has been 
mentioned before that the matter content of the universe behaves like dust during this 
period, leading to a situation in which density fiuctuations grow and could form primordial 
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The length of time for which the massive scalar field dominates is limited by the decay of 
the field into radiation. The Tolman-Bondi models can be used to predict how much time 
is necessary for the collapse to a black hole. In an Einstein-de Sitter universe of density 
Pc = SH^/Sn, a region of dust with density p = pc{l + C) collapses to a singularity in a time 
tc ~ ^ttH^^^^^^'^, where C, is the volume average of (. If the WKB approximation is valid 
then a black hole will form only if the scalar field has not decayed during the time tc- 



The results obtained here for the vahdity of the WKB approximation have been hmited 
to asymptotically fiat spacetimes and time symmetric initial data but they can be adapted to 
remove these limitations for applications to cosmology. We expect there to be little change 
to the results if the length scale A is smaller than the horizon size 1/ H. 

We have a simple condition for validity of the WKB approximation in the Oppenheimer- 
Snyder limit, M^ ^ 1. In the early universe, the mass within the scale A is given by 
M = ^H^\^. The scalar field will behave lilce dust provided that H^X^ji ^ 1. If we tak;e 
H\ ~ 1, then the condition becomes A/x ^ 1, which means that the WKB approximation 
is valid up to the formation of the event horizon if it is valid initially. 
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FIGURES 



Event horizon 




FIG. 1. A section of the Tolman-Bondi metric with coordinates (R, t) for an exponential density 
distribution. The shell r = X and the event horizon are shown, as well as the line t = t* where the 
WKB approximation breaks down.. 




Event horizon 



FIG. 2. A section of the Oppenheimer-Snyder metric with coordinates {R,t). The edge of the 
matter distribution r = A and the event horizon are shown, as well as the line t = t* where the 
WKB approximation breaks down. 
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Critical scalar field mass 
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FIG. 3. Configurations with parameters above the hnes will form black holes. The dotted line 
denotes the top-hat field profile and the solid line denotes the exponential field profile discussed 
in the text, n is the scalar field mass, M the ADM mass and A the length scale. 
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